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Mobility causes network structures to change. In PSNs where underlying network structure is changing rapidly, we 
are interested in studying how information dissemination can be enhanced in a sparse disconnected network where 
.nodes lack the global knowledge about the network. We use beamforming to study the enhancement in the information 
"dissemination process. In order to identify potential beamformers and nodes to which beams should be directed we use 
.the concept of stability. We first predict the stability of a node in the dynamic network using truncated levy walk nature 
'of jump lengths of human mobility and then use this measure to identify beamforming nodes and the nodes to which the 
.beams are directed. We also develop our algorithm such that it does not require any global knowledge about the network 
'and works in a distributed manner. We also show the effect of various parameters such as number of sources, number 
.of packets, mobility parameters, antenna parameters, type of stability used and density of the network on information 
'dissemination in the network. We validate our findings with three validation model, no beamforming, beamforming 
.using different stability measure and when no stability measure is associated but same number of node beamform and 
'the selection of the beamforming nodes is random. Our simulation results show that information dissemination can be 
.enhanced using our algorithm over other models. 

'Keywords: Dynamic Network, Human Mobility, Stability, Information Dissemination, Beamforming 
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"1. Introduction 

In dynamic communication networks such as 
iMANETs, connectivity of the network has always 
'been the point of the discussion because it affects the dis- 
, semination of the information within the network. Recent 
'studies have shown that mobility can both speed up as 
iwell as slow down the dissemination of the information 
'in the network The speed up of the information 

dissemination process is because the nodes have a higher 
'probability of meeting other nodes in the network that 
,either have or do not have the information. While a 
'reduction in the rate of information dissemination process 
could happen as the mobility could leave the network 
disconnected. For example, consider a network with 
one large giant component, the maximum time taken to 
disseminate the information in the network in static case 
is of the order of the diameter of the network. However 
when the network is mobile, the giant component may 
break into smaller components. This disconnection would 
stop the transmission of the information from one com- 
ponent to another until the connection is made or a node 
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with the information connects to the component that 
does not have the information. Other than mobility there 
are many other factors that also affect the information 
dissemination in the network. These include, bursty data 
0], strength of the tie 0], network structure activity 
pattern H, content type @, node characteristics 0, 
altruism [8'] etc. 

Contrary to MANETs, in a static network, 0, 
showed that beamforming at the transmission end could 
make a disconnected wireless network connected and could 
help in faster information dissemination. Though the focus 
of their technique was to achieve small world like charac- 
teristics, their approach also used natural techniques like 
flocking and lateral inhibition to solve the issue of con- 
nectivity in a distributed way without the global knowl- 
edge of the network. Small world properties of a net- 
work are marked by small Average Path Length and high 
Clustering Coefficient {APL was defined as the mean of 
distance between all pairs of nodes in the network while 
CC was defined as the fraction of number of triangles that 
exists to all possible triangles that could exist). Further, 
authors of [y, [l^l used closeness centrality [ll[ to clas- 
sify nodes in the network into 3 categories, 1) nodes that 
have low closeness centrality beamform, 2) nodes to whom 
beams are directed have high closeness centrality and 3) 
nodes that remain omnidirectional. However, in dynamic 
network, it is argued in fl2!| that beamforming can also 
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increase the dissemination of the information. The au- 



thors of showed that only small fraction of nodes are 
needed to beamform for enhancing the information dis- 
semination in the mobile network. They also showed that 
information dissemination process depends on the width 
of the beam and better performance can be achieved if 
beams are longer and narrower if beams are created us- 
ing same power as that of omnidirectional beam. They 
further showed that further gains are achieved if antenna 
rotates. However, in their algorithm they only considered 
single source single packet information dissemination sce- 
nario. Another work related to beamforming in M AN ETs 
include The focus of [13|] was to study the affects of 
beamforming in human centric vehicular networks. The 
authors in [13| also showed that beamforming can enhance 
information dissemination in mobile networks with similar 



and it does not require any global information. In terms of 
degree, if a node is p% stable it would mean that the node 
has maintained p% of its initial set of neighborhood. In 



observations as of [l^l- However the only difference be- 
tween fl2\ and fi3\ being in the selection of nodes that 
beamform. Moreover, contrast to [H used direc- 

tional antenna at reception. 

Moreover, we think that information dissemination can 
can be further enhanced as compared to 12, 1^ in terms 
of time in a dynamic network where connectivity is very 
less. The above statement forms our motivation for this 
paper where we discuss a mechanism to enhance informa- 
tion dissemination in dynamic networks. Further in this 
paper, we consider Pocket S'witched iVetworks, PSNs, a 
type of MANET, where devices are smalL have limited 
energy and are carried by humai is [l4 , 1^ 16 1 . Going back 
to the concept mentioned in llC|. in dynamic networks 
when global data is missing, computation of the central- 
ity measures to determine nodes that should beamform 
in order to reduce APL in a distributed manner is highly 
complex and energy wasteful. Thus the al gori thm in (9i.dQ] 
cannot be used in dynamic conditions. In [17| authors pro- 
posed a mechanism to predict centrality measures. How- 
ever, their method used global data to perform prediction. 
Their method aggregates all the previous snapshot of the 
graphs to perform prediction. Many studies have argued 
that aggregation is not good as it results into false esti- 
mations. Towards this, we provide an efficient way to en- 
hance information dissemination by associating a stability 
measure with a local centrality measure that is computed 
based only upon one previous snapshot and not aggregated 
snapshot. 

A stability measure reveals how stable_a system is. A 
system can be a network 1^ IS IS 21 > 
network 1231 or a node in the network 



22 



a link in the 



Section [121 



provides a brief overview of these methods. As we are in- 
terested in computing node characteristics in the dynamic 
network, stability of network and stability of the link in- 
terests us less. For node stability, a node is said to be 
stable either when it is not moving or when the node is 
moving but it is keeping some percentage of its interac- 
tions intact. In this paper, we would limit our focus to 
the second definition of node stability. Thus we use degree 
to compute stability of a node. Degree is a local measure 



24l | authors used degree over two consecutive snapshots to 



compute stability of the node. In a directional and mobile 
scenario however, computation of the new neighborhood, 
then beamforming and determining again the new neigh- 
borhood would consume energy and time. Thus it would 
be beneficial for a node to estimate its new neighborhood 
and its stability beforehand. This requires a node to be 
able to predict its stability measure at each time instant 
even before knowing about the neighborhood. Towards 
this, the nodes calculate the probability of being connected 
to one of its old neighbor using the mobility characteristics 
of PSN (jump length distribution of human mobility pat- 
tern, i.e., truncated levy walk model). Once the probabil- 
ities are calculated the node average out the probabilities 
to estimate its stability. 

In dynamic network where information dissemination is 
the main issue, low stability of the node would mean that 
the node has met many different nodes. This would in- 
crease the probability for a node to receive the information 
in the network. Moreover, if that node has high degree it 
would imply that the node would be able to receive more 
information through its neighbors and in turn will be able 
to send out more information in the network. Further, ir- 
respective of the degree, if a node is highly stable it would 
mean that the node virtually has same neighborhood and 
therefore has low chances of receiving different informa- 
tion that is flowing in the network. Coming back to our 
original problem of enhancing dissemination of informa- 
tion in the dynamic network, a low stability node becomes 
a prime candidate for a node that has to beamform while 
high stability nodes candidates for nodes toward which 
beams should be directed. This would allow high stabil- 
ity nodes to also receive information in the network. Fur- 
ther, all different scenarios to which nodes can be potential 
beamformers and recipient of the beams are mentioned in 
section [3] 

Thus, in this paper, we use above mentioned discussion 
to build our algorithm and focus on how information dis- 
semination can be further enhanced in the dynamic spatial 
network where nodes lack global knowledge about the net- 
work. We design our algorithm such that it first addresses 
how nodes can predict their stability each time using pre- 
vious instant information and then apply beamforming. 
We then use this stability measure to identify nodes that 
beamform and the nodes towards which the beams should 
be directed in order to disseminate the information in the 
network much faster. Towards this, in section [3] we pro- 
vide an analytical model for predicting the stability of a 
node and provide a mechanism for the nodes to determine 
who will beamform to whom. Extensive simulations led 
us to formulate different scenarios where we study the im- 
pact of various parameters like different number of sources 
and packets, density of the network, mobility parame- 
ters, radius of communication, stability parameters and 
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antenna parameters. We find that information dissemina- 
tion greatly depends on the above mentioned parameters 
and have an impact in percolation of number of nodes hav- 
ing all the packets. This percolation shifts towards left or 
right depending on the parameter used. 

Further, this paper is structured such that it gives a brief 
overview of related work and useful concepts in section [2j 
the model outline in section [3] and the results obtained in 
section m We conclude the paper by presenting conclusion 
of the study, section [S] 



Notation 


Meaning 


a 


power law exponent constant 




exponential cutoff constant 


pO 


probability distribution function 




radius of gyration 


gain 




radiation intensity 


r 


transmission radius 


M 


max. antenna elements available with v 


m 


number of antenna elements used by v to 
beamiorm |m G [2,MJ 


Bb 


bore-sight direction 


G 
K: 


initial network of V vertices and E edges 
a node e V 


t 


time 


K 


neighborhood of Vi at time t 


Gt 


snapshot of G at time t 




adjacency matrix of Gt 


K 


ith eigen value of X* 


t 

Pi 


ith eigen value of X 


vt 


link status between Vi and Vj at t 


S{G) 


stability of the graph 


T 


max. time 


n 


number of words 


I 


max. length of word 


Ent(Vi,Vj) 
Ent{Gt) 


I C 1 ' 1 1 I TT" ITT" 

entropy of a link between Vi and Vj 


entropy of the Graph Gt 


sm 

Aft 


stability of the node Vi 


in-degree neighborhood of Vi at t 


A 


simulation area 


Smax 


max. stability threshold 




min. stability threshold 


h 


jump of neighborhood node of Vi 
jump of Vi 


h 


distance between initial position of 
neighborhood node and the final 
position of Vi 




guarantied connectivity angle 


BWv, 


beam width of node Vi 


BLv, 


beam length of node Vi 


H{x) 


set of high stability nodes in the 
direction x 


P 


density of nodes 



Table 1: Notations and their meaning. 



2. Related Work and Useful Concepts 

In this section, we present related work and some useful 
concepts. We give a brief overview of information dis- 
semination and human mobility model and beamforming 
in section 12.11 12.21 and 12.31 respectively. We then provide 
overview of stability measures in section 12.41 We further 
formulate tabled] that details all the notations used in the 
sections hereafter. 

2.1. Information Dissemination 

Models related to dissemination of information in the 
network have been well studied. These models are primar- 
ily inspired from disease spreading within the population. 
The models mainly use combinations of 4 different states, 
Susceptible, Infected, i?ecovered and _Exposed. A node 
is said to be in state S if it is not infected by the disease 
but is likely to gain the disease in future. State / is when 
the node has the disease. State R corresponds to recovered 
state of the node, meaning that a node was infected before 
but has now gained immunity while E state states that a 
node is exposed to the disease and is highly likely to gain 
infection. One of the main study focusing on the informa- 
tion dissemination in the network was that of Watts and 
Strogatz's in 1998 [25,]. The study infiuenced researchers to 
investigate information dissemination in various different 
types of networks, for example, scale free networks 261.127 1 
and sparse networks [28j]. However, recently the focus has 
been shifted towards the study of information dissemina- 
tion in dynamic networks d [ii iM m m, m, s, . 

In dynamic scenarios, dynamicity affects the dissemi- 
nation of the information in the network [l[. It is thus 
critical while modeling information dissemination to know 
what application is targeted. For example, dissemination 
of a warning signal in the community should be very fast 
however dissemination of viruses should be very slow. In- 
formation dissemination is also affected by many other fac- 
tors like, bursty data Q, strength of the ties Q, network 
structure activity pattern content type Q, node 
characteristics 0, altruism Q, etc. In this paper, we do 
not concentrate on these factors and focus only on mo- 
bility. As we consider PSN, in the next subsection, (Cf. 
section 12. 2p , we would provide some insights to human 
mobility characteristics. 

2.2. Human Mobility Model 

Many mobility models have been proposed that capture 
the characteristics of human mobility. A comprehensive 
survey of mobility properties associated with humans can 
be found in [H, |33| , while a clear difference between dif- 
ferent human mobility models can also be found in [s^. 

Most of the traditional models on human mobility use 
spatial properties of human mobility in the model. How- 
ever, it was brought up that human mobility also has 
temporal characteristics and follows truncated power law 
distribution. Truncated Power Law is characterized by a 
probability distribution that has the properties of a power 
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law with an exponential cutoff, eq. [T] where proportion- 



ality constant is equal to 



0'- 



r(i- 



a is the power law 



decay exponent constant while /3 is the cutoff value. Trun- 
cated power law distribution means that the distribution 
starts as a power law and ends as an exponential curve. 
The exponential cutoff is because 67% of the time human 
movement is limited to a distance called radius of gyra- 
tion, Tg 38|. Other factors that influences this behavior 
are the boundaries assigned to the human movements and 
sampling of the mobility data 0. Studies have revealed 
that a for human movement mainly lies between 1.75 ± 
0.15 while /3 varies [sl]. Figure 1(a) shows the distribu- 
tion of jump lengths of human movements using a = 1.6 



and different /3. While fig. 1(b) shows the distribution of 
jump lengths of human movements using different a but 
/3 = 300.0. The figure 1(a) show that as /? increases the 
jump length also increases. This is because the cutoff for 
large /3's is also high so the jump lengths are high. Further, 



figure 1(b) show that as a increases jump length decreases. 



This is because the distribution of jump length is inversely 
related to a. Further, according to this definition of trun- 
cated power law, the jump length, x, should be greater 
than or equal tol, (a::>l). Aa;<l would increase the 
probability to more than 1. Thus, for a; < 1, p{x) is set to 
1. 



p{x) cx X "e (3 



(1) 



Other properties of human walk that obey truncated 
power law distribution include pause time |40|, inter con- 
tact time t4l| and r^, [sl]. The values of a and /3 for pause 
time are 0.8 ± 0.1 and 17hr respectively. However, these 
values are 1.65 ± 0.15 and 350km for Vg. We next provide 
a brief overview of antenna models and beamforming in 
section [ 



2. 3. Beamforming 

When omni-directional radiations from different 
isotropic and equidistant antenna elements of a node 
interfere in a constructive and destructive manner with 
each other a characteristic long-range beam is achieved. 
The ability of a node to produce such long-range beams is 
called beamforming. Many different types of beamforming 
techniques exist in literature. Some of them are random 
beamforming, greedy beamforming and preferential 
beamforming. A brief survey of beamforming techniques 
used in wireless networks can be found in [IQ]. 

When only single antenna element is used there is no 
interference and thus omni-directional radiation is main- 
tained. However, when multiple antenna elements are 
used, long-range beam are obtained. The gain pattern ob- 
tained through the interference of isotropic and equidistant 
antenna elements is given as eq. [2] as defined by ^] 
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(2) 



where 9 is angle with the z-axis, 4> with the xy-plane 
and u{6, (j)) is the radiation intensity. Research has iden- 
tified many different antenna models using the above 
model. Among them Uniform Linear Array antenna 
model iULA)^ and Uniform Circular Array antenna model 
{UCA) are most common. ULA model is characterized by 
linear deployment of antenna elements while UCA model 
is characterized by circular deployment of antenna ele- 
ments. The gain pattern for the ULA antenna model is 
only dependent on the number of antenna elements and 
has no dependency on the bore-sight direction {B^, the di- 
rection of maximum radiation intensity). Fig. [2] shows the 
gain pattern obtained for ULA model using m = 8 and 
r = 50 where m G M is the number of antenna elements 
used out of M available antenna elements and r is the 
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transmission radius. On the other hand, for the UCA an- 
tenna model, it was shown that gain pattern is dependent 
on both the number of antenna elements and Bb [43] . 



90° 90° 




270° 270° 



(c) = ±60° (d) Bt = ±90°. 

Figure 2: Gain pattern obtained for different Bi,, m = 8 and r = 50 
in the ULA model 



Stability 


Method 


Local/ 


for 




Global 


Graph 


Similarity 


Global 


Graph 


Adj. spectrum [18] 


Global 


Graph 


Adj. matrix [19] 


Global 


Graph 


Adj. matrix [20] 


Global 


Graph 


Rank Overlap [21, 22] 


Global 


Link 


Entropy [23] 


Local 


Node 


Neighborhood [24] 


Local 



Table 2: Related work based on Stability measures. 



edges. We explain more about different stability measures 
in detail in the following subsections. Table [5] shows the 
differences between proposed methods. 

2.4- 1- Similarity 

Similarity is well defined concept in data mining con- 
text. Representing a system as a vector, many similarity 
measures have been defined to determine the similarity 
between two systems. These methods include Euclidian 
Distance, Cosine Similarity, Jaccard Coefficient, Pearson 
correlation, etc [i^, EHj. A high similarity value means 
that two systems are almost similar and no change has oc- 
curred in the system over time. Considering a graph as a 
system, similarity among two instances of a graph can be 
computed easily. As an example. Cosine Similarity is cal- 

^ E ^ n 

culated as — ' , where Ny, is the neighborhood 

of Vi G V and \V\ is the number of nodes in the network 
at time i*. 



We next provide an overview of various stability mea- 
sures available in literature. 

2.4- Stability Measure 

Stability is the measure that determines how much sta- 
ble a system is in a dynamic environment. Stability can 
be calculated by determining the percentage change in the 
system. Stability can be defined at different scales, for 
example, in dynamic graph, stability can be defined for a 
graph, for a node and for a link. Many different models 
have been proposed that can compute stability. Like, in 
[18,] authors use spectral difference to compute the differ- 
ence between two graphs. In authors use adjacency 
matrix to compute the stability of the graph. In [2l|, |22| 
authors rank the nodes based on degree to compute stabil- 
ity of the graph, in |23| authors use entropy to infer stabil- 
ity of the link while in j24| authors define a way to compute 
stability of a node. Moreover, concepts from data mining 
can also be applied to determine the stability as suggested 
by Zayani et al. in [23| . Using data mining concepts, sta- 
bility can be computed using similarity methods between 
two instances of the system into consideration. Consid- 
ering a graphical system, G{V,E), with V nodes and E 



2.4-2. Spectral Distance 

Change in the graphical structure can also be computed 
using the difference in their adjacency spectrum Ad- 
jacency spectrum is defined by the Eigen values of the 
adjacency matrix, X. In [181] . authors use distance be- 
tween adjacency spectrum to compute difference in the 
graph structure. Further, in [18| the authors also survey 
different methods to compute spectral distance. We use 
eq. [3] as defined in 46, I^l to show difference between two 
graphs. 



SD{Gt,Gt') - <^ 



1=1 



\v\ \v\ , 

i=l 1=1 



otherwise 



(3) 

where A* is the i*'' Eigen value of the adjacency matrix of 
the graph at time t while /i* is the ith Eigen value of the 
adjacency matrix of the graph at time t . 
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Stable [191 


Stable [20.1 


1 


1 


Yes 


Yes 


1 





No 


No 





1 


No 


No 








Yes 


No 



Table 3: Truth Table for link between Vi,V,. 



2.4-3. Adjacency Matrix Based Measures 

Stability of a graph can also be measured using changes 
in the adjacency matrix, Hanneke et al. [19| 

proposed to calculate stability of the adjacency as eq. IH 
Considering two possibilities for a link existence, all pos- 
sible status of a link at two given time instances can be 
constructed as table El where X is the adjacency matrix 
and Xy, y. is the status of the link between Vi and Vj at 
t. Hanneke et al. proposed that a link is stable if it had 
occurred in both the instances or did not occur while a 
link is not stable if it changes its state from 1 to or vice 
versa. However, another study of Tang et al. 0| said that 
link with status (0, 0) in Table [3] should not be considered 
for the stability as the link does not exist in any of the 
considered time instances. They proposed that stability 
should be computed as eq. [5] 



SiG) 



J2 ^V,.V,^V,.V, + (1 ~ ^V'.,\/,)(l ~ ^V„Vi) 



\V\{\V\^1) 
SiG) = ^ E ^v-. 



(4) 
(5) 



where 



T 



T-l 



^ yVjev/v, 



V^,VJ 



(6) 



2.4.4- Rank Overlap 

Stability can also be measured by computing rank over- 
lap as suggested in 2l|, |22|. However, done for email net- 
work, the authors suggested to compute degree of all nodes 
in the network at a given time instant. Once the degrees 
are computed, the nodes are ranked based on their degree. 
As the email network was large, top ly = 1000 nodes were 
selected from the ranking. This was done for different 
time instances and lists were made. Commonalities be- 
tween two selections were then found in the form of same 
nodes appearing in the two lists. Mathematically overlap 
was model as eq. [7] as shown in 47 1 



' "i li'iS I \ '!(!'>/ ,1,1 , 



0.8 




100 150 
Timed) 



Figure 3: Relation between different stability measure values. Inset 
shows a zoomed-in image of stability measure computed through eq. 
|4] The stability value using eq. |4]are high and ~ 0.995 because all 
links that do not occur in two consecutive time instances are also 
considered. 



S{G) 



\Gtop n 



Top\ 



(7) 



where G^q^ and G^^p are the list of all top v nodes at 

time t and t respectively. 

Further, we also provide a comparison of stability mea- 
sures as in section [2.4.2112.4.^ and [274 . 41 for graphs. Figure 
[3j shows differences in the stability values for a graph us- 
ing different methods. The stability measures in figure [3] 
are calculated using two consecutive time instance graphs. 
Moreover, rank overlap is computed using top 50% nodes. 
The figure shows that the stability measure differ from 
each other. 

2.4.5. Entropy 

Entropy means the "predictability" of the system. This 
means that Entropy reflects the measure of stability. If the 
system is lacking predictability the states of the system 
will not be well known and the entropy would be high. 
If the entropy is low it would mean that the system is 
predictable. 

Zayani et al. (23j used this definition to predict graph- 
ical systems. Zayani et al. used the entropy with Katz 
centrality to predict the occurrence of the link in the dy- 
namic network. They argued that a link between nodes Vi 
and Vj can be represented as a series of I's and O's collected 
over time, with 1 meaning the existence of the link at a 
time instant and meaning the absence of the link. Za- 
yani et al's algorithm built these sequences only for those 
links which have occurred over time at least once. Their 
algorithm used learning to know about the nodes existing 
in the network. To compute entropy, they used an algo- 
rithm that identified maximal number of unique words, n, 
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Tvs e. z=o 



^00 500 1000 1500 2000 2500 3000 3500 4000 

T 




Figure 4: Variation of £ with respect to T when Z=0 for worst case 
sequence. 



Figure 5: Variation of £ with respect to T when Z=0 for best case 
sequence. 



in the string of I's and O's. They defined entropy of the 
hnk between a pair of node, Vi and Vj, as eq. [5] 



Ent{V„V,) = 



n * hi(rt) 
T 



(8) 



where T is the time scale over which the hnk is ob- 
served or the maximum time. However, Zayani et al 
did not provide any analysis to their work. We per- 
formed a worst and best case analysis of the algorithm 
proposed by Zayani et al. We first identify the worst 
case sequence. Intuitively worst case sequence should 
be "0101010101010101010101....". But, because of the 
periodic nature this sequence has low entropy. Accord- 
ing to the algorithm proposed by Zayani et al, we find 
"01000110110000010100111001011101110000..." sequence 
to be the worst case sequence. This sequence is formu- 
lated using the fact that one can use two unique words 
of length one (0, 1), four unique words of length two (00, 
01, 10, 11), eight unique words of length three (000, 001, 
010, Oil, 100, 101, 110, 111) and so on. The length of this 
sequence can be represented using eq. IHl 



T = 2^+1(^-1) + 2 + Z 



(9) 



where Z is the remainder of 2'+i{i-i)+2 ^ length 
of largest word. The analysis reveals that £ can be given 
by eq. [10] while n in the series can be given by eq. 1111 



W{Y)+\n2 

ln2 
W(Y)+\n2 

ln2 



if z = o 

1 otherwise 



W(V)4-2 In 2 

2^-nr2 - 2 

^ ln2*2 in2 -2*(ty(y)+31n2)+T*ln2 
W(Y)+2\n2 



(10) 
if z = o 

otherwise 
(11) 



where y = [T/2-l]*ln2, Z e Z* and W(.) is the product 
log function. Fig U shows the variation of £ with respect to 
T when Z = 0. It can be seen that the growth of £ is not 
so evident and it almost take T > 3586 to have £ > 8. For 
T < 3, £ < 0, this is contradictory to the definition of £ as 
£ is the word length which cannot be less than zero. Thus 
to compute entropy T should be greater than or equal to 
3 (T > 3). The best case sequence on the other hand is a 
sequence of all I's or all O's. This will lead to eq. [T^] 



T = 



£{£ + !) 



Z 



(12) 



where Z £ [0, £]. Here n = £ when Z — and n ^ £ 
when Z ^ 0. £ is further defined as eq. [13] 



£ = - * (±V8r - 



1) 



(13) 



Variation of £ with respect to T for eq. [T^ is shown in 
fig [5] As £ cannot be negative we just consider £ > 0. If 
above analysis is used, it will reduce the computation of 
words every time as done in the Zayani et al. algorithm 
and can save lot of node's energy. A similar method to 
that of Zayani et al was proposed in [i^ where authors 
used location preferences instead of links. 

Studies have also revealed that entrogy of a g raph can 
also be defined. Various models, [4S: iSO, Isil, |52, |53|, l54| . 
have been developed to compute entropy of the graph. A 
brief survey of entropy measures could be found in fssl 
[56| . Some of above mentioned models, [H^, use Laplacian 
spectrum of the graph to compute the entropy of the graph 
as eq. [H] 
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\v\ 



i?ni(Gt) -$](A*/2)log(A*/2) 



(14) 



\v\ 



^(A*/2)*(l-A*/2) 



where A* is the i Laplacian Eigen value. 

2.4-6. Neighborhood Based 

All the above mentioned measures, except link entropy, 
for computing stability require global knowledge. A simple 
measure for computing node stability in terms of neighbor- 
hood is defined as eq. [15] 



(15) 



in,Vi I 



where 7V*„ y, is the in-degree neighborhood oi Vi £ V at 
time t. The eq. [15] identifies all those nodes that were 
neighbors of the node at time t and t and computes the 
ratio with the actual number of neighbor nodes at t. A 
similar method was proposed in [24*1 . The authors of 
claim to propose a method for the computation of stability 
of the node but the model lacked clarity in the formulation. 

The above model captures the local computation of sta- 
bility measure for the node but requires two snapshots of 
graphs in time, t and t . This would mean that the node 
would first require the knowledge about its new neighbor- 
hood at time t . It would then compute stability, then 
beamform and then recompute the new neighborhood. 
The node would thus require time in performing above 
four steps. It is highly possible that by that time the node 
changes its position. To reduce computation time it would 
be better if the nodes are able to predict stability before. 
In the next section (Cf. section [3]) we propose a stability 
measure that can be computed using local neighborhood 
and does not require new neighborhood for the calculation 
of stability. 

3. Model 

In this section we provide a method to estimate the sta- 
bility of the node in the dynamic environment and how 
this stability could be used toward enhancing the infor- 
mation dissemination in a PSN. We use truncated levy 
walk characteristics of human mobility to build our model. 
Now, let us first consider a given spatial distribution of 
wireless nodes in area A with transmission radius r. This 
will constitute a spatial network, G{V,E), with V nodes 
and E edges. Let us also assume that nodes are mov- 
ing and follow truncated power law distribution of jump 
lengths. When nodes are moving, the network structure is 
constantly changing and affects the connection setup of the 
nodes with the neighbors. Assuming nodes to have com- 
plete information about the network is not good. Thus, 



let us also assume that nodes lack global information and 
only have local information. Further, let us also assume 
that mobility speed is very much less than that of trans- 
mission speed so that when data is being transmitted the 
node do not change their position. We provide nodes in the 
network with two measures, maximum stability threshold, 
Smax, and minimum stability threshold, Smin and equip 
each node with AI antenna elements. The nodes use single 
antenna element for the omnidirectional beam, however, 
uses m G [2, M] antenna elements when beamforming. We 
consider only directional transmission for now and do not 
consider directional reception. Further, we use similar an- 
tenna model to perform beamforming as used in jlOi] . 

Absence of global information will leave nodes to have 
only certain information about the other nodes. Due to 
mobility, we limit this information to node's initial loca- 
tion, node's final position after the jump and the initial 
position of the initial neighborhood nodes. After nodes 
have moved, a node will have new set of neighbors that 
now will act as the initial neighborhood for predicting the 
stability of the node for the next time instant. 

Following the assumptions, we build an analytical model 
for computation of node stability. Let the distance moved 
by the node Vi & V in T = t — t he equal to bi and the 
distance between a neighbor, Vj £ V at T = t, (Vj), and 

V, at T = / , (v/ ) be k where t < t . 

In order to build the model we first assume that the 
node Vi knows how much distance, Ui, its neighbor Vj has 
moved in T = t ~ t. We will then relax this condition. 
The two nodes, Vi and Vj, will be connected only if the 
node Vj moves such that its movement is within 20i and in 
the direction of movement of Vi as shown in Fig. [B] Using 
law of cosine, 9i can be computed as eq. [TH] 



-1 fa l+lf-r^ 



(16) 



The probability of connection between Vi and Vj therefore 
can be given by, ealTTl 



Pviy,{C\ai,k,bi) = — 



(17) 



where C stands for connection. We now remove the condi- 
tional of knowing . In order to remove the conditional of 
knowing a^, we need the distribution of Oj. Due to the fact 
that jump lengths are distributed according to truncated 
power law, this fact can be used to remove the conditional 
of knowing a^. This results into modification of eq. 1171 to 
eq. [IB] 



li+r 



Pv^.y^{C\h,bi) ^ J fA{a^)*Pv^y^{C\a,,kM) (18) 

\h-r\ 

where fA{ai) is the probability distribution function of at 
given by eq. [T] The lower bound of the integral limit. 
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Figure 6: The probability of having a neighborhood node connected 
in the next time instance is given by the relative size of the dark bold 
line. 



I^i — r|, is absolute value because if the jump of node Vi is 
less than r, bi < r, then the lower limit would be r — li, 
however, if jump length is larger than r, bi > r, the limit 
would be k — r. This fact is captured well by the absolute 
value. Moreover, because the jumps cannot be negative, 
i.e. hi < 0, absolute value helps in addressing this fact. 
As discussed in section jump lengths should be grater 
larger than 1 unit. As the units are not defined, for a jump 
length less than 1 , the scale should be changed accordingly 
such that jump length is larger than 1. Also, here li should 
be larger than 0. If li — then cosine term will not be 
defined and the probability will become 0. This fact is 
rather obvious as Vi can not move to the location of the Vj 
as it thinks that Vj is already occupying the same location. 



2 I j2 



COS 



7(a; 



7(b) 



\g{a)) and 
shows the 
shows the 
7(c) shows 



Considering g{a) — — 2^ 
/(a) = fA^ai) * Pvi,Vj{C\ai,li,bi), figure 
variation of possible values of g{a). Fig. 
variation of 6 for different values of k while fig. 
the variation of /(a) on loglog scale for different values of 
li and r with /3 = 300m and a — 1.6. 

As eq. [IH] refers to the probability of keeping the con- 
nection in the next time instant, this probability directly 
relates to the probability of link occurrence in the next 
time instant. The same probability also relates to how 
much stable a connection is. As there can be N*^ y, in- 
degree neighbors of Vi at time t, the total stability of the 
node should thus be given by eq. [TOl 




10" 



1=0.0, r = 30.0 
;=30.0, r = 30.0 

;=ioo.o, r=30.o 

;=300.0, r = 30.0 
;=500.0, r = 30.0 



300 
a [m] 



500 



(a) Variation of g{a) with respect to / and r 




1=0.0, 1=30.0 
;=30.0, r = 30.0 

;=ioo.o, r=30.o 

;=300.0, r = 30.0 
;=500.0. r = 30.0 



300 
a [m] 



(b) Variation of 9 with respect to I and r 




1 = 0.0, r = 30.0 
i=30.0, r = 30.0 
i = 100.0, r=30.0 
i=300.0, r=30.0 
i=500.0. 7- = 30.0 



10^ 



a [m] 



(c) Variation of /(a) with respect to / and r 

Figure 7: Functional analysis with varying I with /3 = 300 and a = 
1.6. All distances are in meters. For I = 0, no graph is displayed 
because cos~^ is not defined at i = 0. 
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(19) 



The probabihty in the eq. HH] is calculated for a link. 
This probability can be very high meaning the link is 
highly likely to exist, ie., highly stable or can be very low 
meaning the link is unstable. Summing of the probabil- 
ities provides the expected value of the number of nodes 
which will still be connected to the node at the next step. 
Here, a large number of highly stable links make the over- 
all node stability high. For example, consider a node Vi 
with in-degree neighborhood at time t as A^*„ y. = 4. Fur- 
ther, consider 3 link probabilities among 4 calculated as 
very high and one as very very low (tending towards 0). 
As the 3 links are highly likely the overall stability of the 
node Vi at t should be high. This fact is captured well 
through eq. [191 

We use this measure of stability to identify nodes that 
are potential beamformers. If the value from eq. [12] is 
less than Smin, the node switches its beam from omnidi- 
rectional to directional beam. As discussed in the intro- 
duction section, the node with low stability is assumed to 
have met more nodes and is likely to have more informa- 
tion. When such a node beamforms to more stable node, 
it is likely to enhance the information dissemination mech- 
anism in the network. However, when a node, Vi, at t has 

WLy, I = but has \Nf^ y^ | > at i' , the S{V) = 0. This 
would mean that it is highly likely that the node does not 
have any information to transmit. It is thus also beneficial 
for another node with low stability to beamform to such 
nodes. Further, when Vi at t has \N^^ v \ ~ ^ ^^"^ ^^^o 

N*„y_\ = at t , the S{Vi) = 0. This would also mean 
that the node is isolated and might not have any informa- 
tion to transmit. It would be beneficial to beamform to 
this node so that this node is not left without the informa- 
tion. As this node can also be the source of information, to 
transmit its information this node must beamform to any 
node in the network. To distinguish between above cases, 
it is necessary to associate the stability measure with the 
degree of the node. We thus, formulate table [Jj that lists 
all possible combinations of stability and degree measures 
for the nodes that would be potential beamformers and 
the nodes to which these beams are to be directed. 

We next describe how nodes that are beamforming can 
determine the angle, width and the length of the beam. 
We use the same adaptation of antenna configuration in 
our model as used in [3, and use 



(20) 



where BLvi = m * r as described for theoretical antenna 
model (Sector model) in jl3| and BWvi the beam width. 
The equation eq. [^Ojcan provide the width of the beam but 
not the direction of the beam and how to determine the 
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Table 4: Beamforming strategies. 

best m. As in [l3| we also use same approach in deter- 
mining the length and the direction of the beam. When 
a node determines that it is a potential beamformer, it 
randomly chooses m € [2, M]. As omnidirectional beam is 
converted to directional beam, the directional beam has to 
sweep „?,T sectors to cover all the directions. The node 
^ |— l 

similar to [9|, [lOj] sweeps all these sectors to determine the 
best sector. The best sector is then determined by 1) num- 
ber of high stability nodes and 2) their stability value. Eq. 
models the same. 



BDv = max 



\^keH{BDt^,) 



(21) 



where H{BDy ) is the set of all high stability nodes in the 
direction BDy and j is a sector £ [1, -gf^]- Once the 
direction is determined we use Sector model to beamform. 
The high stability nodes that were one hop away when the 
beamforming node was not beamforming will affect our 
algorithm. Directing beams to these nodes will not reduce 
the path length but might reduce the neighborhood. Thus, 
we also condition that the high stability node must be at 
least two hops away when the node is beamforming. If 
the node is not able to find any highly stable two hops 
away nodes within its beam length then it reverts back to 
omnidirectional case. 

Further, beamforming makes the network asymmetric. 
A node will be able to transmit the data but will not be 
able to receive the ACK for the transmission made and 
will not be able to know if there is a high stability node in 
the region of directional beam. As we are just concerned 
about the dissemination of the information in a broadcast 
medium, the problem of asymmetry is of less importance. 
For the second problem, again same as 0, [l^] , we assume 
that a highly stable node, just to let the beamforming node 
know about the connection setup, creates the beam in the 
direction of the beamforming node and then reverts back 
to omnidirectional scenario. 
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4. Simulation and Results 

We use a simulation area oi A = 500mx500m to sim- 
ulate our algorithm. The range of average density, p, of 
nodes per unit area is set to 1x10^^. Initially, each node 
operates in omnidirectional mode using m = 1 antenna el- 
ement with the omnidirectional radius as r = 30to. We set 
the maximum number of antenna elements that the nodes 
are equipped with to M = 6. The separation between 
two antenna elements computed using WiFi frequency, 
/ = 2AGhz. We set a = 1.6, P = 300to, Srmn = lxlO-^ 
Smax = 0.7. As ULA model is more realistic we use ULA 
antenna model to simulate our model. Through our sim- 
ulations, we explore the effect time taken to disseminate 
the information in the network. We use Python to simu- 
late our algorithm and use a confidence interval of 95%. 
We average all the results over 50 initial topologies for a 
T = 100. We validate our results with the scenario that 
assumes global knowledge, random beamforming and no 
beamforming. We simulate different information dissemi- 
nation scenarios using simple Susceptible- infected model 
and its variations. Moreover, we show the difference in 
the time taken to disseminate the information in the net- 
work for multiple source multiple packet for ULA antenna 
model and the sector model and also show the affect of var- 
ious other parameters like, mobility parameters, antenna 
parameters and stability parameters on the information 
dissemination in the network. 

4.1. Results 

We simulate six different scenarios for information dis- 
semination, single source having single packet to transmit, 
single source having multiple packets to transmit, single 
source sending update information, multiple sources with 
each having single packet to transmit, multiple sources 
with each having multiple packets to transmit and multi- 
ple sources with multiple packets where sources can join 
over time. As a validation scheme, we compare our model 
results with three cases, namely, when there is no beam- 
forming, when stability measure as proposed in section 
12.4.61 is used and when same number of nodes as found 
by our model beamform but the selection of the nodes is 
random. 

We first show results obtained for the case when there 
is only one source having only one packet to transmit. In 
this scenario, the packet is generated at the start of the 
simulation in this case. Figure 8(a) shows an improve- 
ment achieved when our model is used over the case when 
there was no beamforming done. The results show that 
in almost 26 time steps all nodes in the network have the 
packet when beamforming is done using our model while 
only 86% nodes have the packet after 26 time steps when 
no beamforming is done. The results show high variations 
because the source node is randomly chosen in the net- 
work. In each run, due to difference in the APL of the 
source node to other nodes, each run results into different 
number of nodes receiving packets. The results also show 



improvements over the other two validation models (96% 
nodes have packets after 26 time steps in both the two 
other validation models). Further, the results also prove 
that higher performance is achieved with less energy con- 
sumed when compared to model in section 12.4.61 

Figure |8(b)| also shows an improvement achieved when 
there is only single source having updated packet. The 
initial packet is generated at the start of the simulation in 
this case while the updated packets are generated until 10 
initial time steps. This is done so that percolation could be 
observed. Here, one packet is generated in one time step 
with the packet generated at 10th time step being the most 
updated packet. Here time take to completely disseminate 
the updated packets is shifted by sa 10 time steps than the 
previous case because of the updated packets being gener- 
ated until 10 time steps. All other observations in number 
of nodes having the packets remain same to that of the 
previous case. Further, in this case there is a hill between 
time step and 10 because the updated packets are not 
generated at all time steps between and 10. When at a 
time step an updated packet is not generated, the network 
assumes that there is no updated packet in the network. 
Under same settings, similar results to that of fig. |8(b)[ 
are obtained for the case when there is only single source 
which has multiple packets to disseminate in the network 
(CF. fig. 



Figure 



8(c)) 



8(d) also shows an improvement achieved when 



there are multiple sources which have single packet to dis- 
seminate in the network. We set number of sources to 
be 40 for all the scenarios dealing with multiple sources. 
Similar to the case of single source, the initial packet is 
generated at the start of the simulation by all the sources. 
The results show that in sa 38 time steps all packets can 
be disseminated in the network when beamforming is used, 
however, only « 80% nodes have received all the packets 
when there was no beamforming (CF. fig 8(d)). Further, 
as there are multiple sources, the number of nodes having 
all the packets has a shift in time. It takes « 11 time steps 
for one node to have all the packets from all the sources 
for the beamforming case while sa 16 time steps when not 
beamforming. 



Figure 8(e) also shows an improvement achieved when 



there are multiple sources with multiple packets to dis- 
seminate in the network. Each source randomly generates 
packets until 30 time steps. This makes each source to have 
different number of packets to be transmitted in the net- 
work. This is done so that percolation could be observed. 
The results show that in sa 65 time steps all packets can 
be disseminated in the network when beamforming is used, 
however, only sa 80% nodes have received all the packets 
when there was no beamforming (CF. fig 8(e)). Further, 



as there are multiple sources having multiple packets, the 
number of nodes having all the packets has a shift in time 
when compared to the previous case. This shift is almost 
by w 16 time steps before dissemination takes off. 

Figure [8(f)] shows an improvement achieved when there 
are multiple sources with multiple packets to disseminate 
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in the network where sources can join as time passes by. 
Sources can join only until 20 time steps so that perco- 
lation of the number of nodes receiving the packets can 
be observed. Each source randomly generates packets un- 
til 30 time steps. This makes each source to have differ- 
ent number of packets to be transmitted in the network 
as discussed above. The results show that in « 56 time 
steps all packets can be disseminated in the network when 
beamforming is used, however, only « 75% nodes have 
received all the packets when there was no beamforming 
(CF. fig 8(c) I. Compared to the results obtained for mul- 
tiple source single packet case, (Cf. fig. [8(d)] ), there IS a 
shift of almost ~ 13 time steps before dissemination takes 
off for the beamforming case. 

We next provide results for varying p, a, j3, r, M, Smin, 
Smax and antenna type for the case when there are mul- 
tiple sources with multiple packets to transmit. In these 
set of simulations we show the effect of the above men- 
tioned parameters on information dissemination when our 
model is used. First we show the affect of varying p on 
the model. As p increases the number of time steps taken 
to disseminate the packet reduces. This is because when 
p is high more nodes are in the giant component and the 



connectivity is high, (Cf. fig 9(a) I. 

For varying a it was observed that when a is high the 
time taken to spread the information in the network de- 
creases, (Cf. fig. |9(b)[ ). For varying /? variations in time 
take for all the nodes to receive the information is shown 
in fig. 9(c) For less /3, when the jump length cutoff is 
less, the nodes are not able to move far. This affects the 
nodes meeting different nodes in the network, thus more 
time is required to spread the information in the network. 
While when /3 is more, that is when jump length cutoff is 
high, nodes are able to move far and meet more nodes and 
spread the information in the network faster. 

Next we show the results obtained when r and M varies. 
Results clearly state that when r is high more nodes are 
connected and the dissemination of packets is faster, (CF. 
fig. |9(d)]). For varying M the results state that when M 
is high the information dissemination is much faster, (CF. 
fig. [9(e)| ). This behavior can be explained as when M is 
high the beam length is also high and thus more nodes can 
be connected. 

Variations in Smin and Smax affects the number of nodes 
beamforming and to whom they beamform. Low value of 
Smin reduces the number of nodes beamforming while a 
high value increases that number. This affects the infor- 
mation dissemination in the network, (CF. fig |9(f)| . How- 
ever, low value of Smax increases the number of nodes 
the nodes can beamform to while a high value decreases 
that number. Moreover, as defined in tabled there can be 
other nodes to which beams can be directed that do not 
have high stability. A high value of Smax only reduces the 
weight of a sector defined by eq. HUby some fraction. As 
one sector that has the best weight has to be chosen, Smax 
should not affect the spreading of the information in the 
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network, (CF. fig 9(g) I 



Table 5: Factors affecting dissemination of packet. The table shows 
the effect on the time taken to disseminate information in the net- 
work completely by the increase in the parameter value. 



Next we show results obtained when different types of 
antennas are used. We used two different types of anten- 
nas, namely, U LA and a theoretical antenna model (Sector 
Model). ULA achieves better results than sector model 
because in ULA there is a secondary lobe with substan- 
tial beam length and width, (CF. fig |9(h)| . This mcreases 
the connectivity of the node as compared to when sector 
model is used. 

These results led to tabulate table [5] that summarizes 
the effect of each parameter, described above, on the in- 
formation dissemination. 

5. Conclusion 

From this paper, we have tried to propose a model that 
could enhance information dissemination in the mobile sce- 
narios. The contribution of this paper is four fold. First, 
while doing a literature survey we proposed an enhance- 
ment to Zayani et al's. algorithm. Second, we proposed a 
probabilistic approach to compute node stability in a mo- 
bile environment. Third, based on the stability measures 
we perform beamforming to enhance information dissemi- 
nation in mobile environment. Fourth, the operations are 
fully distributed and does not require any global knowl- 
edge of the network. 

Our results show an improvement in information dis- 
semination when a stability measure is used to beamform 
over three different cases including the case when no beam- 
forming was done and when randomly nodes were chosen 
to beamform. Through our results we have also tried to 
prove that probabilistic calculation of stability can very 
well enhance the information dissemination in the network 
and have validated the results. The probabilistic calcula- 
tion of stability does not require past snapshots but one 
previous one. Further, our algorithm works in fully de- 
centralized manner. The results show that information 
dissemination is affected by number of sources, number 
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(a) % nodes having packets when only one node acts as a source (b) % nodes having packets when only one node acts as a source 

and it has only single packet to transmit. and it has updated packet to transmit. 
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(c) % nodes having packets when single node acts as a source and (d) % nodes having packets when multiple nodes act as source 

it has multiple packets to transmit. and each of them has single packet to transmit. 




Figure 8: Percentage of nodes having packets. 
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(a) % of nodes having packets for variable p. 



(b) % of nodes having packets for variable a. 
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(c) % of nodes having packets for variable /3. 
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(g) % of nodes having packets for variable Smax- (h) % of nodes having packets for different antennas. 

Figure 9: Percentage of nodes having packets when multiple nodej £jct as source and each of them have multiple packets to transmit. 



of packets, mobility parameters, density, radius of com- 
munication, number of antennas, antenna type, stability 
parameters and stability type used. However, in the cur- 
rent work only fixed number of nodes is considered in the 
simulation area. In real scenarios, nodes can join as well 
as can move out of the simulation area. We would like to 
address these issues and see how information dissemina- 
tion is affected due to the dynamic population. Further, 
we would also like to use real mobility traces to validate 
our results. 
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